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Abstract: We calculate the lowest order cross-section for Dirac magnetic monopole produc- 
tion from photon fusion (77) in pp collisions at y/s = 1.96 TeV, pp collisions at \fs = 14 TeV, 
and we compare 77 with Drell-Yan (DY) production. We find the total 77 cross-section is 
comparable with DY at ^/s = 1.96 TeV and dominates DY by a factor > 50 at = 14 TeV. 
We conclude that both the 77 and DY processes allow for a monopole mass limit m > 370 GeV 
based upon the null results of the recent monopole search at the Collider Detector at Fermilab 
(CDF). We also conclude that 77 production is the leading mechanism to be considered for 
direct monopole searches at the Large Hadron Collider (LHC). 



1. Introduction 



Magnetic monopoles have been a theoretical curiosity since the founding of electromagnetic 
theory and have motivated numerous innovative experimental searches. Maxwell's equations 
possess a dual, electric-magnetic symmetry that goes unrealized without the discovery of 
magnetic charges. Stronger motivation for monopole searches was provided by Dirac who 
showed that the existence of a single magnetic monopole is sufficient to explain the observed 
quantization of electric charge Q, an empirical fact which goes otherwise unexplained. Mag- 
netic monopoles are also present in a majority of grand unified theories (GUT) of particle 
interactions where monopoles are produced as topological defects during the GUT phase tran- 
sition 0| . Consequently, if a GUT were realized in the early universe, after an era of inflation, 
then a population of magnetic monopoles would be left over as relics of the Big Bang 0. 

Decades of interest in monopoles has inspired monopole searches in a wide range of 
physical settings ^J. It has been proposed that relativistic monopoles could be observed as 
a component of the cosmic rays |j] and recent flux limits have been reported Q employing 
these techniques. Monopoles searches have been conducted in exotic materials like moon 
rocks [|7j] and terrestrial materials exposed to excessive radiation ||. Collider searches for 
directly produced monopoles have been performed, most recently (9|, |lQ| . Despite all efforts 
to date there are no definitive signals for the existence of magnetic monopoles. 

The quantization of angular momentum for magnetic and electric poles yields the Dirac 
quantization condition, eg = n/2 (setting h = c = 1) for electric and magnetic charges e and 
g, respectively. The magnetic charge 

9 = ^ = (■£-) - (68.5)e (1.1) 
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is large in units of the electric charge (where we choose n = 1 and define a = e 2 ~ 1/137). 
The large monopole charge implies a strong monopole-photon coupling, a characteristic of 
magnetic monopoles that will be exploited in this article. 

A useful theory of monopole interactions does not currently exist to perform direct pro- 
duction calculations. The large monopole-photon coupling precludes the use of perturbation 
theory leaving us with a lowest order approximation as our only means to proceed. Previous 
authors [||, ||, 1C , [ll| have employed a minimal model of monopole interactions which assumes 



a monopole-photon coupling that is proportional to the monopole's induced electric field g(3 
for a monopole moving with velocity (3 = v/c. We follow these authors and use this same 
minimal model in our photon fusion (77) and Drell-Yan (DY) calculations that follow. 

Monopole searches at colliders are restricted to Dirac-type monopoles (and antimonopoles) 
which are hypothesized to be fundamental particles dual to the electron (and positron). 
GUT monopoles are excluded in collider searches as they generally have too large a mass 
(M ~ IOOAgut where Aqut is the GUT symmetry breaking scale) and their internal struc- 
ture exponentially suppresses their production cross-section. Models of Dirac monopole pro- 
duction have relied extensively on the DY process in which a quark and antiquark (qq) from 
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interacting protons annihilate to produce a monopole-antimonopole pair (mm). (For an 
extensive review of the Drell-Yan process see |l2|| .) The CDF Collaboration at the Fermi- 
lab Tevatron recently reported the results of a search for the direct production of magnetic 



monopoles [10|. With no monopole events found, CDF sets a mass limit m > 360 GeV as- 
suming DY production of monopoles. A particle collider that probes a new energy frontier, 
such as the Large Hadron Collider (LHC), will open up new physics possibilities including the 
potential for the discovery of magnetic monopoles. Future monopole searches are likely to be 
undertaken at the LHC and will require detailed simulations of monopole events based upon 
our best knowledge of the leading monopole production mechanisms. Anticipating this need 
we investigate DY production at the LHC and consider an alternative production mechanism, 
the 77 fusion process. 

The 77 production cross-section of heavy leptons, 77 — * L~L + , has been studied in 



comparison with DY in pp collisions at LHC energies |13[. The full 77 process includes the 
individual regimes of inelastic, semi-elastic, and elastic scattering. The lepton production 
cross-section in each regime was found to be of the same order of magnitude while the total 
77 cross-section, the sum of the individual regimes, was found to be nearly 10 2 below the DY 
cross-section. We repeat the 77 calculations for monopole production which entails replacing 
e (for leptons) with gj3 (for monopoles) that, in light of eq. ( |1.1[ ), will lead to greatly enhanced 
cross-sections. The DY cross-section will also be enhanced for monopole relative to lepton 
production, however it remains to be seen which process will dominate. 

We restrict our calculations to Dirac monopoles (hereafter "monopoles"), assumed to be 
spin 1/2 fermions of minimal charge, only consider electromagnetic interactions, and assume 
a monopole-photon coupling g(5 for final state monopoles of velocity j3. In Section (|2|) we 
describe the 77 and DY calculations and discuss their relative dominance for lepton versus 
monopole production. The details of the 77 and DY calculations, which have been widely 
reported in the literature, are presented in Appendix ([A]) for the case of monopole production. 
In Section (H|) we present our results for pp collisions at y/s = 1.96 TeV, discuss the monopole 
lower mass bound reported by CDF, give a mass limit based on our calculations, and present 
our results for pp collisions at y/s = 14 TeV. Our concluding remarks are given in Section 
(!)• 



2. Photon Fusion Versus The Drell— Yan Process for Monopole Production 

The 77 and DY processes for lepton production have been widely studied in proton colli- 
sions. The 77 subprocess for monopole production, depicted in Fig. ([!]), yields a monopole- 
antimonopole pair mm in the final state. The incident photons are radiated from the electric 
charge distribution of the colliding protons (or antiprotons) . During elastic scattering the 
photon couples to the whole proton charge e and during inelastic scattering couples to consti- 
tuant quarks of charge e q = r\e where r\ = 2/3(— 1/3) for q = u,c,t(d,s,b). We will compare 
77 with the DY process (Fig. (|2|)) which dominates for lepton production. The full 77 and 
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Figure 1: The Feynman diagrams for the 77 fusion subproccss which produce a monopole- 
antimonopole pair (mm) in the final state. Two incoming virtual photons (7) are radiated from 
the interacting protons or antiprotons (not shown). The virtual photons couple to the total charge 
distribution of the proton (during elastic scattering, which leaves the proton intact) or to a constituent 
quark within the proton (during inelastic scattering). The monopole-photon coupling is found from 
minimal assumptions of monopole interactions described in the text. The 77 cross-section formula 



for spin 1/2 monopoles is given in eq. (A..3) 




Figure 2: The Feynman diagram for the Drell-Yan subprocess which produces a monopole- 
antimonopole pair (mm) in the final state. An incoming quark and antiquark (qq), constituants 
of the colliding protons, annihilate into a virtual photon (7) which then pair produces the final state 



monopoles. The DY cross-section formula for spin 1/2 monopoles is given in eq. ( A.1C ) 



DY cross-section formulae are presented in Appendix ([A|) with the relevant couplings for 
monopole production. 

To understand the relative strengths of 77 and DY production it is instructive to compare 
the electromagnetic couplings in the total cross-sections. For the estimates to follow we 
only consider quark-photon couplings e q = rje where r] = 2/3( — 1/3) for q = u,c,t(d, s,b). 
The full 77 calculation includes couplings to the whole proton charge which will provide a 
marginal increase to the 77/DY enhancement we find in eq. ( |2.3| ) below. In the case of lepton 
production, 77 suppression relative to DY is anticipated merely by counting the powers of 
electromagnetic couplings in the total cross-sections. The ratio of electromagnetic couplings 
in the lepton production cross-sections for 77 relative to DY is 

4 4 

n = \^ = ffa 2 (2.1) 



where 77 is the average fractional quark charge contributing to the cross-section. Going from 
lepton to monopole production we simply replace e — > g(3 = e(5/2a in the final state couplings 
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(as shown in Figs. (||) and @). In this case the ratio of couplings is 



rm= ) 2 =fy 2.2 

(where a = e 2 ). We can now estimate the change in the 77/DY cross-section ratio expected 
for monopole versus lepton production by taking a ratio of the ratios 

R = = fP/A ^ 47oq 

setting (3 = 1. Drees, et ah, find 77 production of leptons to be nearly 10 2 below DY 
for pp collisions at ^/s = 14 TeV, which implies a factor ~ 50 dominance of 77 over DY for 
monopole production assuming = 1 in eq. (|2.3| ). The effect of /3 < 1 for the production of 
slow moving monopoles is to be determined in our full calculation which follows. 



3. Monopole Production in Proton Collisions 

We calculate 77 fusion for monopoles production following the formalism of Drees, et al. 
[O]. The detailed formulae are presented in Appendix (|A|) and full documentation of our 
77 calculations is reported in a thesis of one of the authors [O]. Unlike the DY process, 77 
production yields equivalent results in pp and pp scattering. The full 77 calculation includes 
contributions from three individual regimes; inelastic, semi-elastic, and elastic scattering, and 
we sum these individual regimes to find the total 77 cross-section. For inelastic scattering, 
pp — ► YY77 — > XX mm, both intermediate photons are radiated from partons (quarks or 
antiquarks) in the colliding protons. To approximate the quark distribution within the proton 



we use the Cteq6-1L parton distribution functions [15] and choose Q 2 = s/4 throughout. 
Following fill] , we employ an equivalent-photon approximation jn| for the photon spectrum of 
the intermediate quarks. In semi-elastic scattering, pp — ► pXjj — > pXmm, one intermediate 
photon is radiated from a quark, as in the inelastic process, while the second photon is 
radiated from the other proton, coupling to the total proton charge and leaving a final state 
proton intact. The photon spectrum associated with the interacting proton must be altered 
from the equivalent-photon approximation for quarks to account for the proton structure. 
To accommodate the proton structure we use the modified equivalent -photon approximation 



of [17]. For elastic scattering, pp — > pp'j-y — ► ppmm, both intermediate photons are radiated 
from the interacting protons leaving both protons intact in the final state. 

The Fermilab Tevatron is &pp collider at center-of-mass energy ^/s = 1.96 TeV. We have 
calculated 77 and DY production of monopoles at the Tevatron and our results are presented 
in Fig. (||). The individual 77 scattering regimes are all of similar magnitude and happen to 
fall near the DY curve. To avoid cluttering Fig. @ the individual regimes are not shown, 
but their contributions to the total 77 cross-section are roughly: inelastic (10%), semi-elastic 
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Figure 3: The total cross-sections for 77 and DY production of monopoles in pp scattering vs. the 
mass of the produced monopole at y/s — 1.96 TeV. The CDF reported cross-section limit at 95% 
confidence level (shown above) can be used to set a lower monopole mass limit assuming a production 
mechanism. The 77 and DY curves are nearly equal at the exclusion limit and independently call for 
a mass limits m > 345 GeV and m > 350 GeV, respectively. We also plot the sum of 77 and DY and 
find a monopole mass limit m > 370 GeV assuming both production mechanisms. 



(50%), and elastic (40%). The total 77 and DY cross-sections shown in Fig. (J3j) happen to 
be nearly equal over the range 300 to 500 GeV with 77 dominating at lower masses. 

The result of a search for the direct production of monopoles was recently reported by 
the CDF Collaboration [10|. The search uses 35.7 pb" 1 of CDF run II data where a specia 
monopole trigger was employed. Monopole event simulations and conservative estimates of 
their experimental acceptance were used to establish a cross-section limit of approximately 
200 femptobarns over a mass range 200 to 700 GeV at 95% confidence level based upon 
a lack of observed monopole events (see Fig. (§)). Assuming DY production of monopoles, 
CDF establishes a monopole mass limit of m > 360 GeV. Their monopole acceptance depends 
upon the production kinematics, but they estimate a limit in the total variation in acceptance 
to be less than 10% and conclude that mass limits from production mechanisms other than 
DY can be set with reasonable accuracy. Thus, we are justified in considering mass limits 
from 77 production based upon the CDF 95%CL limit. 

The DY curve shown in Fig. (3) of crosses the 95% CL limit near 360 GeV while our 
DY curve is slightly lower and crosses near 350 GeV. Therefore, our DY calculation calls for 
a slightly lower mass limit than CDF reports, but the addition of the 77 contribution to the 
DY production argues for an increase in the mass limit of 20 GeV to m > 370 GeV. See the 
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Figure 4: The cross-sections for 77 fusion and DY processes in a pp interaction vs. the mass of 
the produced monopole at y/s — 14 TeV. The 77 process is calculated in the elastic, semi-elastic and 
inelastic regimes shown above and the sum of the individual processes gives the total 77 production 
cross-section (solid line). 



77+DY curve in Fig. @. 

The LHC is designed to produce pp collisions copiously at \fs = 14 TeV. The results of 
our calculations at LHC are presented in Fig. (||). We find that each of the individual 77 
scattering regimes dominates DY by a factor > 10 and the total 77 cross-section is a factor 
> 50 larger than DY. Based upon our results we conclude that 77 fusion will be the leading 
mechanism for direct monopole production at LHC and argue for further investigation of the 
77 process in detailed simulations of LHC monopole events. 

If the LHC were to attain 100 fb _1 of integrated luminosity our calculations predict 
greater than 700,000 monopole events from 77 fusion for a monopole mass of 1 TeV. By 
comparison, the yield of 1 TeV monopoles from DY production is less than 15,000 events over 
the same period of time. The 77 process will allow LHC to extend their monopole search to 
relatively high masses. After collecting 100 fb _1 of data we predict in excess of 50 monopole 
events at monopole masses approaching 3 TeV. 

4. Conclusions 

Motivated by recent monopole searches at the Fermilab Tevatron and the expectation of 
future monopole searches at the LHC we have calculated monopole production from the 
Drell-Yan and 77 fusion processes. We have compared these processes for both pp collsions 
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at y^i = 1.96 TeV and pp collisons at \fs = 14 TeV. Our calculations are limited to a lowest 
order estimate assuming a monopole-photon coupling proportional to the induced electric 
field of a moving monopole. The monopole-photon coupling is strong for monopole velocities 
(3 ~ 1 which has prohibited our use of a pertubative expansion. 

In pp collisions at the Tevatron we showed that the 77 total cross-section is approximately 
equal to DY in the mass range where the CDF collaboration sets a 95% confidence level cross- 
section limit. Based on our results, DY can be used set a mass limit m > 350 GeV and 77 
can be used independently to set a mass limit m > 345 GeV. When both 77 and DY 
production are considered, the sum of the cross-sections implies the monopole mass limit of 
m > 370 GeV. 

In pp collisions at the LHC we found that 77 fusion is the dominant production mechanism 
for magnetic monopoles by more than a factor 50 over the DY process. The inelastic, semi- 
elastic, and elastic regimes each dominate DY by a factor 10 or greater. We conclude that the 
77 process should be considered the leading production mechanism for monopole searches at 
the LHC and emphasize the need for detailed studies of monpole events using the 77 process. 

A. Appendix 

Our 77 calculations follow the formalism and approximations of Drees et al, We cal- 

culate inelastic, semi-elastic, and elastic processes and assume throughout this report that 
the final state monopoles are Dirac-type of minimal charged (n = 1), spin 1/2 fermions, and 
only consider their electromagnetic couplings. The 77 subprocess must satisfy the kinematic 
constraint s = (k\ + k^) 2 > 4m 2 where k\ and k 2 are the virtual photon four-momenta and 
the final state monopole pair has a total rest mass 2m. We assume an effective photon ap- 
proximation [jl6| to describe the photon spectrum of the interacting quark during inelastic 
scattering. The total cross-section for inelastic scattering is 

a p p eL (s) = ^2 [ dxi [ dx 2 [ dzi f dz 2 e 2 q e 2 q , 

1 J 4m 2 1 's J 4m 2 / ' sx\ J 4m 2 / ' sx\Xi J 4m 2 / ' sx\x-zz\ 
■ fq/p(xi, Q 2 ) fq>/p(x 2 , Q 2 )fy/q{zi) f y / q >{z2) d-^(xiX 2 ZlZ2s) (A.l) 

where m is the monopole mass, e q = r]e where r\ = 2/3(— 1/3) for q = u,c, t(d, s, b), and <r 77 is 
the production subprocess cross-section with the center-of-mass energy y§ = yJx\X2Z\Zis. 
The structure function f q / p is the quark density inside the proton and / 7 / g is the equivalent- 
photon spectrum of a quark. We use the Cteq6-1L parameterization of the parton densities 
and chose the scale Q 2 = 3/4. With 

/ 7 /,W = A/fl'W = £ + Z?) ln(Q max /Q mi n) (A.2) 
where Q max = s/4 - m 2 and Q min = 1 GeV 2 . 
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The final state monopole velocity is (3 = (1 — 4m 2 /s) 1 / 2 for the subprocess center-of-mass 
energy y/S. The 77 — > mm total cross-section is 



17(77 ^ mm) 
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In 



(2-/3 2 ) 



(A.3) 



2/3 1-13 

where g A j3 A = (3 4 /lQa 2 using the Dirac quantization condition. The factor a~ 2 will be can- 
celled by two powers of a from eqs. ( |A.2| ) and ( |A.5| ). 

The semi-elastic cross section for pp — ► mfhpX is given by 



^.semi—el. , 



) = 2 y] / dxi / dzi / dz 2 e 2 /< 

„ JAn?/s J Am? j sx\ J Am? j sx\z\ 



'q/p(xi, Q ) 



• f~t( Z l) fj/pfa) Vj~/(X 1 Z 1 Z 2 S) 

The subprocess energy now is given by y/I = yfsx\z\ z 2 . 

cel. 
'l/P 



(A.4) 



For the elastic photon spectrum f el ; (z) we use an analytic expression from [17] given by 



el. 1 



a 
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and where 
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(l + (l-z) 5 



,4 = 1 + 
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0.71(GeV) 2 



Q 
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min 



(A.5) 
(A.6) 



(s + m 2 )(s - zs + m 2 ) 
(s — rrip)^J (s — zs — m 2 ) 2 — 4m 2 zs 



(A.7) 



At high energies Q 2 ^ is approximately m 2 z 2 /(l — z). 

The purely elastic scattering cross-section where both protons remain intact in the final 
state is 



f dz! f 

J Am 2 Is J Am 



dz 2 f* l /p{z\) ff/ p (z 2 ) <7 77 (s = ziz 2 s). (A.S 



• Am? I s J Am? / zis 

In the DY process the annihilating qq pair must satisfy s = (p% + p 2 ) 2 > 4m 2 , for quark 
four-momenta p\ and p 2 , to produce a final state monopole pair of total rest mass 2m. The 
Drell-Yan cross-section for monopole production is 

■1 



„ J Am 2 - 1 s -' 



' Am? jx\s 



dx 2 f q /p{xi) fq/ p (x 2 ) a q g(s = xix 2 s) 



Iq/p 



for the DY subprocess 



&{qq — > mm) 



12s 



2- V 



(A.c 



(A.10) 



where 77 is the fractional quark charge in units of e and the quark sum ranges from t,b, ...,b,t, 
ensuring that only quarks and antiquarks of the same flavor contribute. 
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